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ABSTRACT 

The recent Fermi detection of the globular cluster (GC) 47 Tucanae high- 
lighted the importance of modeling collective gamma-ray emission of millisecond 
pulsars (MSPs) in GCs. Steady flux from such populations is also expected in the 
very high energy (VHE) domain covered by ground-based Cherenkov telescopes. 
We present pulsed curvature radiation (CR) as well as unpulsed inverse Compton 
(IC) calculations for an ensemble of MSPs in the GCs 47 Tucanae and Terzan 5. 
We demonstrate that the CR from these GCs should be easily detectable for 
Fermi, while constraints on the total number of MSPs and the nebular B-field 
may be derived using the IC flux components. 
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1. Introduction 

CO 

o 

Q\ \ The launch of Fermi Gamma-ray Space Telescope (Fermi) on 11 June 2008 heralded 

a new era for high-energy (HE) Ga mma-ray Astronomy The Large Area Telescope (LAT) 



O 



aboard Fermi (lAtwood et al.l 120091 ) is ~ 25 times more sensitive than its predecessor, the 
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Energetic Gamma Ray Experiment Telescope (EGRET). Recently, Ferm i anounced the de - 
tection of the globular cluster (GC) 47 Tucanae as a bright point source (lAbdo et al.ll2009l ). 
This emission may plausib ly be due to the large number o f millisecond pulsars (MSPs) 



residing in this cluster (e.g. lHarding. Usov. &: Muslimovl 120051 . hereafter HUM05) 



GCs are furthermore possible sources of very high energy (VHE) gamma rays (IBednarek fc Sitarek 



20071 . hereafter BS07), and are / will be important targets for gro und-based Cherenkov t ele- 
scopes such as the High Energy Stereoscopic System (H.E.S.S. - lAharonian et al.l 120061 ) as 
well as future telescopes such as the Cherenkov Telescope Array and the Advanced 

Gamma-ray Imaging System (AGIsjA 



Venter fc de Jager J2008ah modeled the collective curvature radiation (CR) expected 



from 47 Tucanae, demonstrating that this pulsed component should be detected by Fermi, 
depe nding on the numbe r of visi ble gamma-ray pulsars in the GC MSP population, iV vis (see 



Venter fc de Jagerl (j2008bl ) extended these results by modeling the unpulsed inverse 



Compton (IC) and synchrotron radiation (SR) flux components expected due to the interac- 
tion of HE electrons, ejected from MSP magnetospheres, with the ambient cosmic microwave 
background (CMB) and bright starlight photons, in the presence of a nebular magnetic field 
B. 



In this Letter, we give updated CR spectra ($31 see I Venter & de Jagerl (120091 ) for details), 
in addition to further IC calculations (§__J), for 47 Tucanae and also Terzan 5, another GC 
containing many MSPs (§2]). The roles of Fermi, H.E.S.S., and CTA for deriving model 
constraints are discussed in 851 



2. Globular Cluster Millisecond Pulsars 

GCs are very old galactic sub-structures, with evolved stellar composition. As a con- 
sequence, they are likely to harbor a larger than usual density of end products of stellar 
evolution - compact objects. This is supported by observations of X-ray bin ary systems 
and M SP populations in GCs. Since the discovery of the first GC MSP in M28 (ILyne et al. 
19871 ). a total of 140 pulsars have b een discovered in 26 GCs (lFreirdl2009l ). with 90% having 
periods P < 20 ms (IRansoml 120081 ). 



Among the GCs, 47 Tucanae and Terzan 5 stand out by the number of radio-detected 



1 www. ct a-observatory. org 

2 http://gammal. astro. ucla.edu/agis/index.php/Main_Page 
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MSPs they host: 23 and 33 respectively (jRansoml 120081 ) . However, these GCs are other- 
wise remarkable for different reasons. 47 Tucanae is one of the most massive GCs, after 
uj Centauri, while the core densi ty of Terzan 5 gives it the highest expected rate of binary 
interactions (jPooley fc Hutll2006l ). Basic properties for these two objects are summarized in 
Table ffl 

Besides pulsed emission, detected already in radio, X-rays, and gamma rays, MSPs are 
also expected to produce steady low-flux VHE emission. Given the sensitivity of current 
instruments, it is likely that only a population of such objects would produce detectable 
VHE fluxes O- The extension of the cores of the GCs is furthermore clearly below the 
angular resolution of current HE and VHE instruments. In this Letter we are therefore 
not concerned about individual MSP fluxes, but attempt to constrain a particular class of 
polar cap (PC) pulsar models using the total flux from a GC MSP ensemble. This approach 
reduces uncertainties in the CR and IC fluxes due to unknown pulsar geometries, as the 
relative uncertainties in ense mble-averaged fluxes, whic h scale as the reciprocal of the square 
root of the number of MSPs (IVenter fc de Jagerll2008al ). are typically one order of magnitude 
smaller (for ~ 100 MSPs) than for the case of single MSP 'geometry-averaged' fluxes. 

An important issue in the context of this work is the uncertainty on the size of the 
MSP population in each system, beyond the radio-detected ones. We dinstinguish between 
two numbers. We define the number of visible gamma-ray pulsars (iV v ; s ) as those MSPs 
whose CR is (mostly) beamed toward earth, so that observers will in principle be able to 
measure pulsed emission from these (although it might not always be possible to distinguish 
individual members due to experimental limitations). In contrast, N tot represents the total 
number of MSPs in the GC, whether visible in pulsed gamma rays or not. 

Mostly off-axis pulsars might still contribute to the total CR flux, as observers may be 
clipping the faint edges of gamma-ray beams from these pulsars. We therefore include them 
in A/vis, along with the on-beam MSPs. HUM05 noted that gamma-ray beams of MSPs are 
probably wider than their radio beams, so iV v i S is probably larger when applied to gamma-ray 
pulsars (as done here) than when applied to radio ones, possibly up to N vis ~ N tot . The IC 
flux however depends on N tot , as all pulsars, including ones with off-beam geometries, are 
expected to contribute VHE leptons (mostly electrons, as most MSPs are believed to be pair 
production starved; see HUM05) which upscatter soft photons to gamma-ray energies. 

The isolated location of 47 Tucanae in the So uthern sky away f r om th e Galactic Plane 
made it a prime target for GC studies. In radio, iMcConnell et al.l (120041 ) estimated from 
the unresolved emission a maximum of 30 detectable MSPs. In X-ray, Chandra observations 



revea led a large number of unidentified sources in the core of 47 Tucanae (IHeinke et al. 



20051 ) . from which followed an upper limit of about 60 MSPs, based on the X-ray properties 
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of radio-confirmed MSPs in the GC. 



For Terzan 5 , the available observations are not so constraining. Again in the radio, 
Fruchter &: Gosd (120001 ) estimated a range of 60 — 200 MSPs in the core of the GC, from Very 
Large Array (VLA) and Australia Telescope Compact Array (ATCA) data. No equ i valent 
result has yet been derived in X-rays. The analysis of Chandra data by iHeinke et al.l (120061 ) 
yielded X-ray sources, but with too high a sensitivity limit to constrain the MSP population. 



Finally, for both systems, GC evolution models (e.g. Ilvanova. Heinke. fc Rasid 120081 ) 
suggest MSP numbers below 200, of which a fraction would be detectable. 



3. Pulsed Gamma-ray Flux 



As previously (jVenter fc de Jage r1 l2008ah . we use an isolated pulsar PC model (IVenter fc de 



1997 



2005) , including the effect of Gene ral Relativistic (GR) frame dragging (e.g. lMuslimov fc Harding 



model ( 


Venter & de 


fager 


Muslimov & Harding 





Harding &: Muslimovlll998l ). In accordance with the expectation that most MSPs are 
inefficiently screened by CR and IC scattering pairs (HUM05), we use an unscreened ac- 
celeration potenti al. We furthermore us e the same population of MSPs: P-values (period- 
derivatives) from iBogdanov et al.l (120061 ) . which imply an average spin-down luminosity of 
E TOt = 2.2 x 10 34 ergs -1 , as well as canonical values for equation of state (EOS) parameters 
(moment of inertia Ins — OAMR 2 , stellar radius R = 10 6 cm, and pulsar mass M = 1.4M Q ). 
As previously, we add the CR spectra of N vis = 100 MSPs, randomly selected from our 
population of 13 members (each with a random pulsar geometry, i.e. magnetic inclination 
and observer angle, including off-beam cases) iteratively a million times, and so obtain an 
average pulsed spectrum, with 2<r-bands indicating the uncertainty due to varying P and 
P-values, as well as the unknown pulsar geometries. 



Venter fc de Jagerl (j2008al ) used a delta function approximation for the CR spectrum 
radiated per prima ry, but the resulting; spectr al shape did not correspond well to results of 



similar studies (e.g. lFrackowiak &: Rudakl|2005l ). as pointed out by Dyks (2008, personal com- 



munication). We therefore developed a corrective technique involving the full CR spectrum 
per primary (I Venter fc de Jagerl l2009i ). and applied it here. 



Assuming that the MSPs in Terzan 5 have similar basic pulsar parameters than those in 
47 Tucanae, we scaled the cumulative CR spectrum of the latter to distances of d = 5.5 kpc 
and 8.7 kpc. The resulting updated CR spectra for 47 Tucanae and Terzan 5 are shown in 
panels (a) and (b) of Figure [TJ 
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4. Unpulsed Gamma-ray Flux 

We calculate the ejection spectrum of electrons leaving each MSP by binning the num- 
ber of primary electrons ejected per unit time according to their residual energies at the 
light cylinder. We add the spectra of N tot = 100 MSPs (with random inclination angles) 
and iterate the Monte Carlo procedure a million ti mes. From t his we obtain an average 



cumulative ejection particle spectrum and 2<r-bands (j Venter! 120081 ). 



We divide the region where unpulsed radiation is generated into two zones: Z0, reaching 
from r = to r = r c , with r c the core radius, and Zl, reaching from r = r c to r = r^ m , 
with rhm the half mass radius (see Tabled]). For 47 Tucanae, we used energy densities of 
^rad ~ 3 000 eV/cm 3 for Z0, and u] &A ~ 100 eV/cm 3 for Zl for the st arlight component, 



assum ing; a temperature of T = 4 500 K (more details may be found in IVenter &: de Jager 



2008bl ). These become u® ad ~ 1000 eV/cm 3 and u] &d ~ 40 eV/cm 3 for Terzan 5 (note that 
the energy densities, and therefore the IC 'bumps' associated with the upscattering of stellar 
photons, are linearly dependent on the assumed total stellar luminosity - see Tabled]). For 
the CMB component, we use mcmb ~ 0.27 eV/cm 3 (for T = 2.76 K). 

We next calculate the steady-state injection electron spectrum for each zone by multi- 
plying the cumulative ejection spectrum by an effective time-scale which incorporates both 
radiation losses and particle escape. We modeled the latter (i.e. trapping of particles by the 
nebular field) assuming Bohm diffusion. The injection spectra were calculated for nebular 
fields of B — 1/iG and B = 10/xG (BS07). Using these results, we obtained the expected 
unpulsed IC flux for each GC (see Figured]). While the lower-energy SR spectra are not dis- 
cussed here, we do include SR losses in the calculation of the steady-state injection particle 
spectrum (see IVenter Sz de Jager l2008bh . 



5. Discussion and Conclusions 

We have modeled the HE and VHE components of CR and IC flux expected from a 
population of MSPs in the GCs 47 Tucanae and Terzan 5, for 100 members (Figured]). For 
illustration, we assume that all MSPs are visible (i.e. N vis = N tot = 100). However, the CR 
spectra are linearly dependent on iY" vis , while the IC spectra are linearly dependent on iV tot . 



Fermi has recently released a list of bright sources (> lOa. lAbdo et al.ll2009l ). including 
preliminary fluxes of 47 Tucanae of (5.1 ±2.0) x 10~ 8 cm _2 s~ 1 and (5.6 ±1.0) x 10~ 9 cm _2 s _1 
for energy bands spanning 100 MeV - 1 GeV and 1 GeV - 100 GeV respectively. Our average 
CR spectrum is overpredicting these fluxes by factors of ~ 2 and ~ 7. However, the pulsed 
spectrum may quite reasonably be scaled to A^ vis = 50 (as suggested by observations). Also, 
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extrapolating from observed trends in single MSP spectra, one may speculate that it is 
possible that the cut-off energy may be smaller due to uncertainties in EOS parameters, 
the electric potential, magnetospheric structure, and pulsar geometry. For instance, taking 
iV vis = 50, and halving all CR spectral energies (i.e. shifting the CR spectrum to the left by 
a factor of 2), the model predictions improve to ~ 0.5 and ~ 1.7 times the Fermi flux in 
these bands. 

Our ensemble- averaged CR spectrum, resulting from 10 s single MSP spectra (i.e. 10 6 
100-member cumulative spectra) with typical spectral cut-offs around a few GeV (generally 
depending on the individual MSP's P, P, EOS, and geometry), exhibits a broad peak around 
~ 3 GeV. This average spectral cut-off energy follows directly from the assumed acceleration 
potential and B-field geometry of the PC model under consideration, as well as the fact 
that we simulate a random selection from the P — P parameter space for GC MSPs (similar 
to HUM05), and with different geometries. Given the fact that our CR spectrum seems 
to be reasonably successful at reproducing the preliminary Fermi data, we expect that the 
predicted cut-off energy must be close to the real one. This implies that the average model 
acceleration potential must be of the right order of magnitude, assuming dipolar B-fields. A 
more detailed study involving further model and parameter space constraints will have to 
await further Fermi spectral results. 



We indicate EGRE T upper limits (ULs) from lFierro et al.l (119951 ) on Figure la, and from 



Michelson et al.l (119941 ) on Figure lb, converted assuming an E~ 2 spectrum. Interestingly, 
the 100 MeV UL is close to the Fermi sensitivity for 47 Tucanae. The EGRET UL for 
Terzan 5 is not very constraining, and implies iV v i S < 363 for d = 5.5 kpc, and iV v i S < 907 
for d = 8.7 kpc (using the average integral CR spectrum). 

As a reference, we have also included a CR prediction from HUM05, scaled to N vis = 100, 
in Figured) HUM05 noted that the collective radiation from MSPs in GCs may be visible for 
Fermi, and proceeded to apply their unscreened GR frame- dragging PC model to the MSPs 
in 47 Tucanae. They estimated the CR flux above 100 MeV from 47 Tucanae by integrating 
along a single magnetospheric B-field line, normalising to a (conservative) solid angle of 1 sr 
per PC, and scaling by the standard Goldreich- Julian PC current. This approach, which was 
followed to circumvent full 3D modeling, may be regarded as yielding a theoretical upper 
limit, given that different pulsar geometries (especially off-beam ones) have not expicitly 
been taken into account as done here, thus significantly exceeding our results at 100 MeV. 

Predictions from BS07 focused on the possibility that HE leptons, which may be ac- 
celerated inside MSP magnetospheres and / or reaccelerated at colliding pulsar and stellar 
wind shocks, gradually diffuse through the GC and upscatter CMB and soft stellar photons 
to produce an unpulsed component in the GeV-TeV energy band. BS07 assumed that 1% 



- 7- 



of the spin-down power (with the population average taken as E TOt ~ 1.2 x 10 35 ergs" 1 ) is 
injected in relativistic leptons. In the absence of detailed model predictions, they assumed 
power-law injection spectra (typical for shock acceleration) with various low and high-energy 
cut-offs and spectral indices, and normalized these using the assumed average MSP spin- 
down power and conversion efficiency Representative resulting IC spectra, linearly scaled 
to our population-average spin-down luminosity of E rot = 2.2 x 10 34 ergs -1 , are shown in 
Figure [TJ Our IC predictions, obtained by conservatively assuming no reacceleration of elec- 
trons, roughly agree with those of BS07 around a few TeV, but deviate at other energies due 
to different injection spectra, as well as slightly different assumptions regarding the energy 
density of stellar photons in the GCs (compare our two- zone approach with BS07's energy 
density profile given in their Figure 1). 

Constraints on N tot as well as B may be derived from our model, using the predicted 
IC flux and assuming non-detection above 1 TeV at the level of 1% and 0.1% of the Crab 
flux, as achievable by H.E.S.S. and CTA, respectively. As shown in Figure [2j regions of the 
B — N t ot parameter space could be excluded from VHE ULs, e.g. choosing B = 10/iG for 
47 Tucanae, the H.E.S.S. UL constrains the total number of MSPs to < 33 (and < 3 for 
CTA). Terzan 5 is less constrained (these limits become < 141 and < 14 for d = 5.5 kpc, 
and < 353 and < 35 at d = 8.7 kpc). For small B, particle escape is greater (and IC flux is 
lower), so that the constraints are less severe. This is also true for large B, where higher SR 
losses lead to inhibited IC flux. Furthermore, taking reacceleration of VHE electrons into 
account would increase SR and IC fluxes, implying more stringent constraints on iV tot and 
B in the case of a non-detection. 

The B — N^t constraints are however dependent on the choice of diffusion coefficient 
(a larger coefficient would result in enhanced particle escape and thus lower IC flux). It 
is thus possible to invert the argument and obtain constraints on the diffusion coefficient 
in the absence of reacceleration, using future VHE limits and assuming reasonable values 
for iVtot and B. To this end, a combination of Fermi data and Monte Carlo modeling 
may help constrain N Y \ R and N v \ R /N tCf t., and thu s N to t- Also, measurements of the predicted 



diffuse ultraviolet SR (jVenter fc de Jagerll2008bl ) and unpulsed diffuse gamma-ray IC spectra 



may further constrain the lepton injection spectrum as well as B. Lastly, if B is large 
enough (> 20 /iG) so that diffusion losses are unimportant relative to SR losses, particle 
reacceleration in the GC may be studied without having the concern of the latter being 
masked by energy- dependent diffusion. 



We would like to thank Alice Harding, Dave Tompson, and Jarek Dyks for useful discus- 
sions. This research was supported by an appointment to the NASA Postdoctoral Program 
at the Goddard Space Flight Center, administered by Oak Ridge Associated Universities 



- 8- 



through a contract with NASA, and also by the South African National Research Founda- 
tion and the SA Centre for High Performance Computing. 



REFERENCES 

Abdo, A. A. et al. (for the Fermi-LAT Collaboration), 2009, submitted to ApJS 
(|arXiv:0902.1340p 

Aharonian, F. et al. (for the H.E.S.S. Collaboration), 2006, A&A, 457, 899 

Atwood, W. B. (for the Fermi-LAT Collaboration), 2009, submitted to ApJ 



(arXiv:0902.1089) 



Bednarek, W., & Sitarek, J. 2007, MNRAS, 256 (BS07) 

Bogdanov, S., Grindlay, J. E., Heinke, C. O., Camilo, F., Freire, P. C. C, & Becker, W., 
2006, ApJ, 646, 1104 



Castellani, M., 2009, A galactic globular clusters database, http://snipurl.com/gclusters 
accessed 19 Feb. 2009 

Cohn, H. N., Lugger, P. M., Grindlay, J. E., & Edmonds, P. D., 2002, ApJ, 571, 818 

Frackowiak, M., & Rudak, B., 2005, Adv. Space Res., 35, 1152 

Fierro, J. M. et al, 1995, ApJ, 447, 807 

Freire, P. C, 2009, Pulsars in globular clusters, www.naic.edu/~pfreire/GCpsr.html, ac- 
cessed 19 Feb. 2009 

Fruchter, A. S., & Goss, W. M., 2000, ApJ, 536, 865 

Harding, A. K., Usov, V. V., & Muslimov, A. G. 2005, ApJ, 622, 531 (HUM05) 
Harding, A. K., & Muslimov, A. G. 1998, ApJ, 508, 328 
Harris, W.E., 1996, AJ, 112, 1487 



Heinke, C. O., Grindlay, J. E., Edmonds, P. D., Cohn, H. N., Lugger, P. M., Camilo, F., 
Bogdanov, S., Freire, P. C, 2005, ApJ, 625, 796 

Heinke, CO., Wijnands, R., Cohn, H. N., Lugger, P. M., Grindlay, J. E., Pooley, D., & 
Lewin, W. H. G., 2006, ApJ, 651, 1098 



- 9- 



Hinton, J. A., 2004, New Astron. Rev., 48, 331 

Ivanova, N., Heinke, C. O., & Rasio, F. A., 2008, in 40 Years of Pulsars: Millisecond Pulsars, 
Magnetars and More, ed. C. Bassa, Z. Wang, A. Cumming, & V. M. Kaspi, AIP Conf. 
Ser., 983, 442 

Lyne, A. G., Brinklow, A., Middleditch, J., Kulkarni, S. R., Backer, D. C, & Clifton, T. R., 
1987, Nature, 328, 399 

McConnell, D., Deshpande, A. A., Connors, T., & Abies, J. C, 2004, MNRAS, 348, 1409 
McLaughlin, D. E., & Van der Marel, R. P., 2005, ApJS, 161, 304 

McLaughlin, D. E., Anderson, J., Meylan, C, Gebhardt, K., Pryor, C, Minniti, D., & 
Phinney, S., 2006, ApJS, 166, 249 

Michelson, P. F. et al., 1994, ApJ, 435, 218 

Muslimov, A. C, & Harding, A. K. 1997, ApJ, 485, 735 

Ortolani, S., & Barbuy, B., Bica, E., Zoccali, M., & Renzini, A., 2007, A&A, 470, 1043 
Pooley, D., & Hut, P., 2006, ApJ, 646, L143 

Ransom, S. M. 2008, in 40 Years of Pulsars: Millisecond Pulsars, Magnetars and More, ed. 
C. Bassa, Z. Wang, A. Cumming, & V. M. Kaspi, AIP Conf. Ser., AIP Conf. Ser., 
983, 415 

Venter, C, 2008, Ph.D. Thesis, North- West University, Potchefstroom Campus 

Venter, C, & de Jager, O.C., 2005, ApJ, 619, L167 

Venter, C, & de Jager, O.C., 2008a, ApJ, 680, L125 

Venter, C, & de Jager, O.C., 2008b, AIP Conf. Ser., 1085, 277 

Venter, C, & de Jager, O.C., 2009, in preparation 



This preprint was prepared with the AAS IATgX macros v5.2. 



-10- 



Table 1. Selected Globular Cluster Characteristics 





47 Tucanae 


Reference 


Terzan 5 


Reference 


Distance d (kpc) 


4.0 


1 


5.5 - 8.7 


2, 3 


Mass (M ) 


1.0 x 10 6 


4 


3.5 x 10 5 


5 


Core radius r c (arcmin) 


0.40 


6 


0.18 


6 


Half mass radius r^ m (arcmin) 


2.79 


6 


0.83 


6 


Detected radio MSPs 


23 


7 


33 


7 


Total stellar luminosity (L & ) 


7.5 x 10 5 


8 


1.5 x 10 5 


8 



Re ferences. 
(3) ICohn et al- 



ii iMcLaugh 



(12002) : (4) 



in et all (120061): (2) lOrtolani et al. 



Mc 



Ivanoya. Heinke. fc Rasiol (120081 ): (6) 
f jl996h : (7) lFreire! J2009h : (8) BS07. 



(I2007h : 

.aughlin fc Van der Marell (l2005h : (5) 

Castellanil J2009I ). extended from Harris 
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Fig. 1. — Differential CR and IC spectra for 47 Tucanae (la) and Terzan 5 (lb). The 
top IC spectra are for a nebula r field B = 10 nG, w hile the botto m ones are fo r B = 1 fiG. 



Dark gray lines indicate Fermi ( lAtwood et al.ll2009l ) and H.E.S.S. (lHintonll2004l ) sensitivities 
(appropriate for each GC). Arrows represent EGRET upper limits. In (la), we indicate a 
scaled prediction from HUM05, as well as 2cr-bands (light gray bands) due to unknown 
pulsar geometry (the bottom CR one stops where la exceeds the mean value). For both 
panels, we show predictions from BS07 (scaled to an average spin-down luminosity of E rot = 
2.2 x 10 34 ergs -1 ). In (lb), solid lines are for d = 5.5 kpc and dot-dashed ones for d = 8.7 kpc 
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Fig. 2. — Constraints on the total number of MSPs N tot and nebular magnetic field B. 
Panel (2a) is for 47 Tucanae, (2b) for Terzan 5. The dark gray indicate regions excluded 
by a H.E.S.S. non-detection, while light (plus dark) gray regions are excluded by a CTA 
non-detection (see text for details). In panel (2b), the solid lines are for d = 5.5 kpc, and 
dashed lines for d = 8.7 kpc. The thick horizontal lines indicate the number of currently 
detected radio MSPs (see Tabled]). 



